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I. INTRODUCTION 



Various cosmological observations, mainly Cosmic Microwave Background (CMB) jl[ 2], 
Supernova type la (SNIa) 3, 4], Weak Lensing j^j, Baryon Acoustic Oscillations (BAO) 6], 
and redshift surveys such as 2dF Galaxy Redshift Survey (2dFGRS) j3] at low redshift and 
DEEP2 redshift survey 8J at high redshift, have provided cross-checked data to determine 
cosmological parameters with high precision. These parameters imply that our approxi- 
mately 13.72 Gyear-old universe is nearly spatially flat, homogeneous and isotropic at large 



scale, i.e. a Friedmann-Lemaitre-Robertson- Walker (F. 
and has entered an accelerating phase since z rs 0.46 
called Concordance Model (or ACDM model), the universe consists of 4.6% baryonic matter, 
22.8% non-relativistic unknown matter, namely dark matter (DM), and a remarkable amount 
of 72.6% smoothly distributed dominant component, dubbed dark energy (DE) 2]. 

The equation of state (EoS), w = p/p, of DE, is the main parameter which determines 
the gravitational effect of DE on the evolution of the universe, and can be measured from 



1/RW) universe with zero curvature, 
4|. Moreover, according to this so- 



observations without need to have a definite model of DE. Strong evidences imply that the 
EoS of DE lies in a narrow range around w ~ — 1 and has a smooth evolution (4, 9|. However, 
slow variation of w in cosmic time is not excluded. Theoretically, one can classify the EoS of 



101 ] . namely the phantom divide line (PDL). That is, 



DE with respect to the barrier w — — 1 
DE with the EoS of w = — 1 employs the cosmological constant, A, with a constant energy 



density. The case with dynamical EoS where w > — 1, is referred to as quintessence 11] and 



12 



13|, 



w < — 1 corresponds to an odd theoretical matter case, known as phantom energy 
which potentially has elusive properties. For instance, phantom energy density increases with 
time as the universe expands, which consequently leads to divergence of the scale factor in 
finite time 13[. Moreover, sound speed in this kind of matter may be superluminal. In 



addition to these three subclasses, there is another suggestion for the EoS of DE, the so- 
called quintom [ijj], in which w crosses the phantom divide line as time evolves. It behaves 
like quintessence in some periods of history of the universe, and like phantom in the rest. 
It is worth noting that a single scalar field with conventional features minimally coupled 
to gravity cannot cross the barrier according to a corresponding no-go theorem 151 ] . For a 
complete discussion about the origin of DE, see for example Refs. jiol . li| and references 
therein. 
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In addition to the above classification according to EoS, one can consider different models 
for DE. Actually, there is a wide variety of models concerning dynamics for DE; for instance, 
quintessence models which exploit canonical fields with dynamical EoS 111 ], k-essence models 



with a non-canonical kinetic term 
a quintessence and one as a phantom 
dark energy via a superffuid 



17 1 , q uintom models usually with two scalar fields, one as 



14|, Chaplygin gas models which unify dark matter and 



tachyon fields motivated by string theory [19(, chameleonic 
fields in which a scalar field is coupled to matter [20, |2l|, and modified gravity models 
regarded as modifications to Einstein-Hilbert action [22]. 

The modified gravity models are obtained in different ways, for example, by replacing 
Ricci scalar, R, by f(R) as the well-known f(R) modified gravity [23|, by adding additional 



scalar terms like the Gauss-Bonnet term 



24( , or by non-minimally coup 



the Ricci scalar such as the Brans-Dicke (BD) theory of gravitation 25| or its generalization 



ing of a scalar field to 



(GBD) by employing the BD parameter, as a function of cosmic time. 

To explain the late time acceleration of the universe as mentioned it is most often the case 
that the scalar fields interact with matter; directly due to a matter Lagrangian coupling, 
indirectly through a coupling to the Ricci scalar or as the result of quantum loop corrections 
26| 31 ]. If the scalar field self-interactions are negligible, then the experimental bounds on 
such a field are very strong; requiring it to either couple to matter much more weakly than 



gravity does, or to be very heavy 



32|-35|. Unfortunately, such scalar field is usually very 



light and its coupling to matter should be tuned to extremely to small values in order not 
to be conflict with the Equivalence Principal 2fj]. 

An alternative attempt to overcome the problem with light scalar fields has been suggested 
in chameleon cosmology j^llj^]. In the proposed model, a scalar field couples to matter 
with gravitational strength, in harmony with general expectations from string theory whilst 
at the same time remaining very light on cosmological scales. The scalar field which is 
very light on cosmological scales is permitted to couple to matter much more strongly than 
gravity does, and yet still satisfies the current experimental and observational constraints. 
The cosmological value of such a field evolves over Hubble time-scales and could potentially 



cause the late-time acceleration of our Universe 



38 



The crucial feature that these models 



possess are that the mass of the scalar field depends on the local background matter density. 



While the idea of a density-dependent mass term is not new 
in |2l| [3J 



44|, the work presented 



is novel in that the scalar field can couple directly to matter with gravitational 
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strength. 



211 ]. we investigate an 



In this work in addition to the most expectations presented in 
integration of GBD and chameleon models, calling it as chameleonic generalized Brans-Dicke 
model (CGBD), which allows the scalar field to couple with both gravity and matter and 
the BD parameter to be a function of the scalar field. Since the scalar field, which treats 
very lightly on large scales, can couple to matter much more stronger than gravity does, 
the observational constraints can be satisfied by this model. As mentioned, to construct 
a model which its EoS can smoothly cross the PDL, it is necessary to add unconditional 
features to scalar field in order to have extra degrees of freedom. We show that the CGBD 
model studied in this work is capable to predict not only a PDL crossing scenario but also 
a bouncing behavior for the universe. 

To directly measure the expansion history of the universe, it is necessary to find a probe 
which is model-independent and is based on fewer assumptions. This can be performed by 
focusing on dynamics of homogeneous and isotropic large scale features of the universe, that 
is, the Robertson- Walker (RW) geometry. An example of such a probe is the cosmological 
redshift drift (CRD) first proposed by Sandage and McVittie 45|. This method 46| measures 
the drift of redshift of specific absorption lines, specially those arising from the Lyman-a 
transition in the spectra of high-redshift galaxies and quasars (QSOs) in a time interval 



of several years 



47| . The drift, z, maps the Hubble parameter, H(z), and consequently 



determines acceleration or deceleration of the universe, and only vanishes when the expansion 
rate of the scale factor is constant. However, for varying expansion, it is of order of the 
Hubble parameter, approximately 10~ 10 , and gives rise to an acceleration of order of few 
centimeters per second per year, and can be swamped by peculiar acceleration of cosmological 
objects. A second probe is the observations of the difference between the distance modulus 
- red-shift relation for SNIa that verifies the late-time accelerated expansion of the universe 
and will be discussed in here. 

The manuscript is arranged as follows. In the next section, we introduce the CGBD model 
and analytically study conditions for the EoS parameter to pass the PDL. In Section III, 
by numerical calculation a possible bouncing of the universe and EoS parameter crossing is 
examined. In Section IV, the inclusion of observational data in our model through CRD and 
the difference in the Distance Modulus (fi{z)) tests are investigated. Section V represents 
the results and summary. 
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II. THE CGBD MODEL 



We consider the CGBD gravity in the presence of cold dark matter with the action given 



by, 



S = J d^x^f^g ( 



+ M)L r . 



(1) 



where R is the Ricci scalar and u>(4>) is the scalar field dependent BD parameter. The 
matter Lagrangian is modified as f(<f))L m , where /(</>) is an analytic function of the scalar 
field. This term represents the nonminimal interaction between the matter and the scalar 
field. Metric variation of action (JTT) gives, 



G 



[IIS 



1 



f{4>), 



(2) 



In a flat FLRW cosmology, field equations (j2J) become 

o a 2 p m f{4>) Q «0 . w(0)^ 2 



- 3- 



a- 



and 



a a 



a (j) 



(3) 



(4) 



a a* <p a<p 2 <p 2 <p 

On the other hand, the variation of action ([!]) with respect to the spatially homogeneous 



scalar field (p(t) yields, 



a 3 + 2w(0) 



(Pm ~ 3p ro ) ( 7(0) " \<t>f (0)) " ^'(0)^ 



(5) 



where the prime denotes derivative with respect to 0. In the above equations, p m and p m 
are respectively energy density and isotropic pressure of the barotropic perfect fluid. The 
EoS of the baryonic and dark matter is p m = w m p m . From equations ©-(ED , one can easily 
derive the modified conservation equation 



— + 3H(1 + w m ) = -3w m t 

Pm J 



(6) 



which readily integrates to yield 



M 



71 



(7) 
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M 



where M is a constant of integration. The equation ([7]) shows that for / = 1 ,p m 
as expected p m decreases with increasing f(4>). In comparison with the standard Friedmann 
equations one may identify p eS and p cS from equations (J3J) and @ as 



and 



3H 2 



2H + 3iT 



Pmf 0,,„ 



Pmf . 1 - 3w m p m 1 

l 20 3 + 2w0 ; ' 



(8) 



(9) 



where we have supposed p eff = wp cH . Consequently conservation equation reads as 



PcB +3H(l + w)p ca = 0. 



(10) 



From equation ([3]) we obtain, 



H 



± 



'p m / , 3 + 2w< 



30 



One can rewrite the second Friedmann equation, 



V 6 d> 2 <h 



Differentiating ffl2|) with respect to time gives, 



12 



as 



[1 + 3lU m )p m f 



H = — 



3u 



H 



^r^-(^+ir)T-(2^-7T)T 



'l + 3w r 



-Pmf 



20 



(11) 



(12) 



3Wm)(Pm/) 

2(f) 



(13) 



To have the PDL crossing, that isw^Oasw— > — 1, one should have H = —\{l+w)p cfi = 
and = — \wp cS — \{\ + w)p cS ^ Q when w — — 1. From equation (113"]) . we observe that 
it is always possible to have non vanishing and also vanishing H when w = — 1, if one of 
the following conditions is satisfied: 

1) ^ 0, = 0, w m = -1/3, (p m f) = 0, 

2) ^ 0, = 0, iu m = -1/3, (p m f) = 0, 

3) w m -1/3, = 0, = 0, (p m /) = 0, 



4) (p m f) ^ 0, = 0, 



0, w r 



-1/3. 
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In the first case, (1), H(t) is given by equation (ITT]) , and 



H 




t 






(14) 



By imposing the condition H, we obtain, 



H 




uj{1 + 2cu) 







(15) 



Therefore the conditions for having PDL crossing is that i) u ^ 0, ii) u ^ —1/2, or iii) 
Co 2ujH and also the expression inside the bracket does not vanish. In addition, we need 



In the second case, (2), H(t) = ±* H = and H = — ^V. The conditions for 



since one of the PDL crossing conditions is not satisfied, the crossing does not occur. In 
case (4) again since H = 0, the crossing does not occur. 

Finally, the Hubble parameter as a function of redshift, H(z), is obtained using ([8]) 
and (fTOj) as 



which by numerical methods can be used to trace the Hubble parameter versus redshift. 
Evidently, combination of the result of this model-dependent relation with the model-free 
observational data such as CRD and distance modulus, can be used for verification of model. 
With the analytic discussion for the PDL crossing given in this section, a numerical compu- 
tation will be performed in the following section to illustrate the bouncing and PDL crossing 
dynamics for the model. 

III. BOUNCING BEHAVIOR AND PDL CROSSING 

We explore the possible bouncing of the universe. To examine this, choose t — to 
be the bouncing point. Our effective energy density and pressure, equations (IE]), (EJ) and 
consequently derived u shows that it approaches negative values less than —1 in two periods 
of time. Also by numerically solving the set of friedmann equations and scalar field equation 
we obtain H(t) that provides a dynamical universe with contraction for t < 0, bouncing at 



^0. 




H 2 (z) = e^<- 




(16) 
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t = and then expansion for t > 0. The above analysis clearly can be seen in the numerical 
calculation given in Fig. 1. 




t 



Fig. 1: The graphs of the scale factor a and the Hubble parameter H, plotted as functions 
of time for u = tuo<j)(t) n , co = -1000, n = -2, / = f Q e^\ fo = -2000, and b = -1. 
Initial values are 0(0) = —2 , 0(0) = 1 and d(0) = 0.1. 




Fig. 2: The graphs of EoS parameter w, plotted as functions of time for 
u = uj <j){t) n , n = -2, / = foe^, b = -1. For u = -1000 and four different 
values of f (left panel) and for f = —1000 and four different values of u (right panel). 
Initial values are 0(0) = -2, 0(0) = 1, d(0) = 0.1. 

As Fig. 2 shows, in our model, the EoS parameter crosses PDL with w > 0, at different 
times, depending on the initial values and model parametrization. In addition, with respect 



S 



to the bouncing, during the contracting phase, the scale factor a(t) is decreasing, that is, 
a < 0, and in the expanding phase we have a > 0. At the bouncing point, a = 0, and so 
around this point a > for a period of time. Equivalently, in the bouncing cosmology, the 
Hubble parameter H runs across zero from H < to H > 0, and at bouncing point H = 0. 
This shows that a successful bounce requires H > around the bouncing point. According 
to equation (fl2|) . this result is equivalent to 

#=-~(l + ™)p eff >0. (17) 

From Fig. 2, we see that at t — > 0, w < — 1 and H > which satisfy the above condition. 
One of the well known problems in standard cosmology is the big bang singularity where the 
universe with zero volume has infinite density. We see that in our model at the bouncing 
point the scale factor a(t) is not zero and thus we avoid singularity faced in the standard 
cosmology. 

We note that Fig. 2 just shows the dynamics of the EoS parameter as a function of time. 
At early time, crossing from the phantom phase to the non phantom phase towards w = 
looks like a transition say from inflation era to matter dominated era in early universe. The 
model does not show otherwise. For a more observationally oriented perspective, in the 
next section two cosmological tests for observable quantities are performed in terms of the 
redshift. 



IV. THE COSMOLOGICAL TESTS 



A. CRD Test 



As mentioned in section II, the relation ( fl6l) maps the expansion history of the universe 
for a given EoS. On the other hand, the effective EoS lies in a narrow range around w = —1. 
Following 48(, in Chevallier-Polarski-Linder (CPL) parametrization model one can use linear 
approximation, 

z 



dw cpi 

Wq -Ha 

da 



(a - 1) = Wq + U>1 



where w is current value of the effective EoS, and w\ 



l + z- 

dw cp i 



(18) 



is its running factor. Using 
Eq. ( TlBl) we can obtain the following equation for Hubble parameter 



H(z 



^- = On(l + Zf + (1 - n m )(l + Z )3(l+»0+m) x 



H 2 



cxp 



-3wi( 



l + z' 



(19) 



In CPL model the parametrization is fitted for different values of wq, w\ and Q m . In addition, 
in FLRW geometry, the CRD can be extracted from 



[l + z)H -H(z), 



(20) 



where immediately leads to velocity drift 

v = cH 



cH(z) 
l + z 



(21) 



The ve 
Fig. 3 



ocity drift with respect to the source for CPL model against redshift is shown in 



491 ] . A comparison of the model with the observational data shows that the best fit 



values are for wq = —2 and w± = 3.5 in the region z > 3. 



0.4 



-0.2 



V -0-4 



-1.0 



-1.2 




w =-2.5, Wj = 3 



-2.1, Wj = 3.2 



Fig. 3: The graph of the CPL parameterization model in comparison with 
the observational data, plotted as functions of z for different values of wq and w\ 

To draw the left panels of Figs. 4, 5 and 6, we have used the effective EoS parameter 



w for some plausible functions f(4>) and w (</>). The right panel of Fig. 4 shows the vei 



drift against z for different values of fo, with respect to the observational data 



47 



49|. 



ocity 
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z z 

Fig. 4: The graphs of the EoS parameter and velocity drift v, plotted as functions 
of z for tu = UQ(f)(z) n , n = —2, / = foe b ^ z \ and b = — 1 for four different value of /o- 
Initial values are H{0) = 1, 0(0) = -10, 0(0) = 3. 

Alternatively, one may take a fixed fo and changes uq and numerically calculate the EoS 
parameter and velocity drift (see Fig. 5). 




Initial values are H(0) = 1, 0(0) = -10, 0(0) = 3. 
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In another attempt one may consider a case where both fo and ujq change similar to wq 
and w\ in the CPL model. 



w 




V -0.4 




-1.0 



Fig. 6: The graph of the scale factor a and velocity drift v, plotted as function 
of z for to = ujQ(f)(z) n , n = —2, / = foe b ^ and b = —1 for four different values of Uq and fo. 
Initial values are H{0) = 1, 0(0) = -10, 0(0) = 3. 



B. The difference in the distance modulus, fi(z) 

The difference between the absolute and apparent luminosity of a distance object is given 
by, /i(z) = 25 + 5 log 10 di(z) where the Luminosity distance quantity, di{z) is given by 

^ = < 1+ *>jf#j- (22) 

In our model, from numerical computation one can obtain H(z) which can be used to 
evaluate fi(z). To best fit the model for the parameters /o , Uq, H , b, n and the initial 
conditions 0(0), and 0(0) with the most recent observational data, the Type la supernovea 
(SNe la), we employe the y 2 method. We constrain the parameters including the initial 
conditions by minimizing the \ 2 function given as 

XW/O, "o, i, n, m ,m) ^ f WW-^g-M-^)-^)-^' , (23) 

i=l 0i 

where the sum is over the SNe la sample. In relation (1231 . fij he and fi° bs are the distance 
modulus parameters obtained from our model and from observation, respectively, and a 

12 



is the estimated error of the n° bs . In our model the best fit values occur at /o = — 7, 
luo = 1.2, H = 0.84, b = -.4, n = -2, 0(0) = 1.5, 0(0) = 1, with xLn = 804.87232. In 
Fig. 7, the distance modulus, fi(z), in our model is compared with the observational data 
for the obtained parameters and initial conditions using x 2 method. Also in Fig. 8 the two 
dimensional likelihood distribution for model parameters /o , and uq, H$ are shown. 




Fig. 7: The graphs of the distance modulus fi(z) and EoS parameter w 
plotted as function of the redshift for the best fitted model parameters 
f = _7, Uo = 1.2, H = 0.84, b = -.4, n = -2, and ICs. 0(0) = 1.5 0(0) = 1. 
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Fig. 8: The graphs of two dimensional likelihood distribution for parameters fo, Hq and 
u , H with b = -.4, n = -2 and ICs. 0(0) = 1.5, 0(0) = 1 and d L {0) = 0. 



V. RESULTS AND SUMMARY 

In this paper, we have studied the evolution of gravitational and scalar fields in CGBD 
cosmological model in which a light scalar field nonminimally coupled both to the matter 
Lagrangian (as chameleon field) and the scalar curvature (as brans-Dicke field) in the action. 
The solution shows that the evolution of the scale factor of the universe is non-singular in 
a bouncing scenario, with an initial contracting phase which lasts until to a non-vanishing 
minimal radius is reached and then smoothly transits into an expanding phase. Also the 
dynamics of the EoS parameter with respect to time in early universe is shown in Fig. 2 with 
a transition from phantom phase ( say phantom inflation) to non phantom phase towards 
matter dominated era. 

We then analyze the model with the CRD test. The variation of velocity drift for different 
values of fo and ujo with redshift z is shown in Figs. 4, 5 and 6. In addition, the evolution 
of the cosmological EoS parameter in terms of the redshift, with a transition from phantom 
phase (w < — 1) in the past to the non phantom phase (w > — 1) in the recent past is favored 
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for values of fo and uq, in agreement with the current observational data. The PDL crossing 
occurs for different values of fo and ujq within the range of observationally accepted redshift 
z. 

A comparison between our model with the CPL parametrization model shows that our 
model is in better agreement with the data for values of redshift z greater than 3. To best 
fit model parameters to data we need more redshift velocity drift data within the range 
z < 2.5. The PDL crossing in Figs. 4 and 5 shows that in the former the crossing occurs 
for —0.06 < z < 0.03 whereas in the latter for 0.03 < z < 0.06. In regards to the drift 
velocity with respect to the observational data it shows that both velocity drifts are in 
better agreement compare to the CPL model for z > 3 whereas the velocity drift against 
z in Fig. 5 is even more satisfactory. Fig. 4 and 6 show that today and in near future the 
universe experience PDL crossing, in future it will be in phantom phase whereas in near 
past was in non phantom phase. In comparison, Fig. 5 shows that in future and today it 
is in phantom phase, whereas PDL crossing occurred in very near past and in the past the 
universe was in non phantom phase. 

In the second cosmological test, the difference in the distance modulus fi(z) is calculated 
for the model in comparison with the most recent observational data, the Type la supernovea 
(SNe la). To best fit our model with the data we use the x 2 an d likelihood fitting test. The 
X 2 fitting in the model depends on the parameters fo, Uq, H q , b, n and the initial conditions 
(p(0), and 0(0). The multi-dimensional fitting analysis allows us to constrain the parameters 
that result in a minimum x 2 ■ From Fig. 7 left), we observe that both our model best fitted 
distance modulus and ACDM model are compared with the observational data. Although 
both models best fitted with the data for z < 0.5, our model better fitted for larger redshifts. 
In the Fig. 7 right) which the distance modulus is fitted with the most recent observational 
data by using x 2 method, we have a more observationally reliable EoS parameter. The 
graph shows that the current value of the EoS parameter is about ~ 1.1. In future the 
universe approaches phantom phase but later transit to the non phantom phase. In the 
near past the universe EoS parameter crossed the PDL, and approaches the non phantom 
phase in larger redshift. In Fig. 8, the two dimensional likelihood distribution for model 
parameters are given. The graphs show that, the model parameters are bounded in the 
regions — 8 < fo < 05, 0.6 < uj < 1.8 and 0.6 < H < 1.2 with the best fit values are —7, 
1.2 and .84, respectively. 
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In figures 1-6 as can be seen, using different initial values for the model parameters, the 
trajectories eventually converge. So different initial values for the model parameters may 
initially cause divergence among the trajectories but eventually do not effect the dynamics. 
With respect to the initial conditions, comparing the graphs of EoS parameter in Figs. 4, 
5 and 6 (with ICs. 0(0) = —10, 0(0) = 3), with the best fitted EoS parameter in graph 7 
(with ICs. 0(0) = 1.5, 0(0) = 1), one concludes that since the EoS parameters in different 
scenarios show similar behavior, the model is not sensitive to the initial conditions. 

In summary, we may conclude that the CGBD model can predict the PDL crossing and 
fit the observational data for velocity drift and distance modulus better that the CPL and 
ACDM models respectively, subject to constraints on the model parameters. 
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